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^ In this study, we examine the dynamics of a one-dimensional Frenkel-Kontorova chain consisting 

. of nanosize clusters (the "particles") and photochromic molecules (the "bonds"), and being subjected 

to a periodic substrate potential. Whether the whole chain should be running or be locked depends 
^ ', on both the frequency and the wavelength of the light (keeping the other parameters fixed), as 

■ observed through numerical simulation. In the locked state, the particles are bound at the bottom 
of the external potential and vibrate backwards and forwards at a constant amplitude. In the running 
state, the initially fed energy is transformed into directed motion as a whole. It is of interest to 
note that the driving energy is introduced to the system by the irradiation of light, and the driven 
mechanism is based on the dynamical competition between the inherent lengths of the moving object 
(the chain) and the supporting carrier (the isotropic surface). However, the most important is that 

^ ' the light-induced conformational changes of the chromophore lead to the time-and-space dependence 

of the rest lengths of the bonds. 

PACS numbers: 05.45.-a, 66.90. -fr, 63.20.Ry 

The driven dynamics of a system of interacting particles in atomic or mesoscopic scale has been attracting much 
' attention. It is very important in technology, very rich in physics, and widely applicable in many fields, such as mass 
^ [ transport, conductivity, tribology, etc. As far as we know, one of the most typical examples conducted on the driven 
dynamics has been the driven Frenkel-Kontorova-type systems biased by dc, and ac external forces, respectively. 
The Frenkel-Kontorova (FK) model ||l| may describe, for example, a closely packed row of atoms in crystals, a layer 
Q • of atoms adsorbed on crystal surface, a chain of ions in a "channel" of quasi-one-dimensional conductors, hydrogen 
^ ' atoms in hydrogen-bonded systems, and so on. In general, such a system can be treated with two parts in the driven 
dynamical model: the moving object and the supporting carrier. First, the moving object is considered an atomic 
^ ■ subsystem, in which the interparticle interaction is taken as a harmonic interaction between the nearest neighbors. 
^ , Second, the action of the supporting carrier to the moving object is modeled as an external potential, a damping 
constant, and a thermal bath. When, for instance, an external dc driving force / is applied to such a system, its 
^\ response can be very nonlinear and complex. The overdamped case (ry 3> wq, where rj is the external damping and loq 
: is the vibrational frequency at the bottom of the periodic potential) , has been studied in a number of papers . 
When the driving force / changes, the multistep and the hysteretic transition of the system from the locked state to 
the shding state in the underdamped case (77 ^ loq) has also been delineated in detail H,^. At the same time, various 
intermediate regimes can be described by resonance phenomena Q and by the moving quasiparticle excitation, kinks 

^ ®- . 

A theoretical demonstration of driven dynamical behavior, as a steady-state response to an applied driving force, 
must be well verified in laboratory. It is of interest to note that a Josephson junction array is just an technically and 
I I physically excellent demonstration, of many typical nonlinear features of the driven FK-type systems biased by, for 

■ instance, dc forces and ac forces Q. 

C In this paper, we propose a new photo-driven Frenkel-Kontorova model, which may be actualized in experiment. 

^ The basic challenge of our proposal is to keep the time-and-space dependence of the rest lengths of the harmonic 
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oscillator, as was suggested originally by Proto et al. |11| when discussing the atomic scale engines. In this model 



system, the supporting carrier is taken as an isotropic surface, and the moving object (a FK chain) consists of nanosize 
clusters (the "particles") and photochromic molecules |Q (the "bonds"). The particles are identical with mass m; the 
bonds are fiexible with elastic strength k and photoactive with sensitive coefficient 6; the time-and-space dependence 
of the free equilibrium rest lengths is due to the light-induced conformational changes of the chromophore. In contrast 
to the particles, the bond's mass is so small that it can be neglected. For simplicity we restrict the discussion to one 
dimension at zero temperature. The N equations of motion read as 
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where Xn is the coordinate of the nth particle with 1 < n < A^. The second term of Eq.(jl|) describes the dissipative 
interaction (damping) between the particles and the surface. It is proportional to the relative velocities of the particles, 
with proportionality constant 77. The static interaction between the particles and the surface is represented by the 
periodic substrate potential 



$ (x„) = $0 1 - cos 
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Concerning the interparticle interaction, we take the form of the nearest neighboring harmonic interaction 

k 2 
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Due to the light sensitivity, the bond lengths change with the irradiation of light. The free equilibrium rest lengths 
an,ra+5 (i) depend on both the bond's position, specified by the indices n, n + (5, {8 — ±1), and time t. We restrict 
ourselves to a certain choice for an^n+& (t)- 
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where b [b < 1) is an absorption coefficient (amplitude) less than one, oq is the rest length of the harmonic interaction, 
and A2 and tu are the wavelength and frequency of external irradiated light, respectively. Obviously, oq and b are 
intrinsic parameters, and A2 and uj are externally adjustable. Light energy, if pumped into the system in an irradiating 
manner, may produce spatially and temporally correlated changes of the lengths a„.„+5 (t). Then, the dynamical local 
competition between the periodicity Ai and the rest lengths a„.„+5 (t) may induce a driving power, which will impel 
the chain into motion in some manner. Granted, this power must be strong enough to overcome the resistance. 
After introducing dimensionless parameters a;„ — s- (27r/Ai)x„, t (27r/Ai) ^o/mt, and the renormalized param- 
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equation reads as follows: 
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In the present survey, the mass of the particle is to = 1, the period and amplitude of the external sinusoidal potential 
are Ai = 2tt and $0 = 1, respectively. So, if we take the natural unit like this, i.e., the time scaled by (27r/Ai) ^/Wo/m, 
the length scaled by (27r/Ai), and the mass scaled by (1/m), then the physical quantities, time, velocity and so on, 
become dimensionless. And one should, in order to restore them, multiply spatial variables by (Ai/27r), times by 
(Ai/27r) a/to/$o, rnass by to, and energies by $0, etc. 

In our model, parameters 77, fc, 6, A2, w are all very important. For example, when the external viscous damping 
coefficient "q is overdamped rj loq, underdamped rj <ti loq, or other intermediate regions, the driven system exhibits 
very different dynamical behavior, which has been studied in a number of papers as mentioned above. Note that 
in the present survey, the characteristic frequency of particle vibration at the minimum of the external potential is 
t^o = 1- To simplify the problem, we keep the damping coefficient and the elastic strength fixed as 77 — 0.1 and 
k = 1. Here, the value of k also corresponds to the intermediate case between a strong coupled (Sine-Gordon) and a 
weakly coupled chain. Thus, we focus on the parameters b, A2 and lo, all of which definitely characterize how the free 
equilibrium rest lengths an,n+5 (t) depend on the time and position of the bonds. 
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Our simulation of the Runge-Kutta iteration method starts when the particles are at the bottom of the substrate 
potential well at rest, and the boundary condition is chosen to be the simplest, as Xj+^r = Xj+2'KN (periodic boundary 
condition). In this way, ao = Ai = 27r. Only a finite system possible, its size should be appropriate to ensure that M 
waves {M is the wave-number) can be inserted into the background structure of iV-particles system; the integers N 
and M must satisfy the relation: N/M = A2/A1. In all calculations, we took the natural unit (m = 1, Ai = 27r and 
$0 = 1) a-iid kept 77 = 0.1 and fc = 1, as mentioned above. 



In Fig.l we show the velocity of the center of mass of the system Vc [vc = ( X]i=i ) / ( ) — (X^i ^i) /^j 



irii = m] as a function of time t for some different parameters. The results tell us the following: (1) After passing 
through a transitory relaxation, the system turns quickly into a steady-state; (2) There are two types of steady-states, 
the locked state (Fig.l(a)-(c)) and the running state (Fig.l(d)-(e)). In the locked state, the particles are bound at the 
bottom of the substrate potential and vibrate backwards and forwards at a constant amplitude. In the running state, 
the initially fed energy is transformed into directed motion as a whole, with almost constant speed except for a little 
undulation. With different parameters, there are different oscillating amplitude (in the locked state) or running speed 
(in the running state). On the whole, the system exhibits interesting and parameter-dependent dynamical behavior. 

The amplitude b intrinsically characterizes the light sensitivity of the photochromic molecules (bonds) and deter- 
mines the elastic energy stored. To find more about it, we started from 6 = 0, slowly increased it by small steps 
A6 ~ 0.01. With every b value, we let the system relax to 10^ time steps, 10^"^ each, to ensure a stationary state. We 
then calculated the average velocity of the centre of mass (vc) over the time period of 10^ ~ 10^ steps. The results for 
A2/A1 = 2.5, 3 and w = 1, 3, respectively, are shown in Fig. 2 (a) and (b). The processes both have a {vc) increasing 
from zero as soon as the amplitude exceeds a critical value be- At zero speed, each particle is locked in one external 
potential well without any observable macroscopic motion, expecting its emancipation as the driving energy exceeds 
a certain value. This is easy to understand for a dissipation system. 

It is no surprise as well, that each particular situation is unique, given different frequencies and wavelengths as 
temporal and spatial modulations. To get further information of their effects, our next simulation was focused on 
the external photo-driving parameters. Figure 3 presents (vc) as a function of the driving frequencies uj for a fixed 
value of the wavelength A2. From these numerical results, three things are note- worthy. First, as the magnitude of 
frequency slowly increases from w = 0, a dynamical transition occurs, from an initially locked state to a running state, 
and finally back to a locked state again. Second, there exits a peak of maximal velocity. Third, as the wavelength 
increases, the running state has a lower maximal velocity and exists in a narrower span of frequency with lower values. 
The same is true for (vc) as a function of wavelength, as is presented in Fig. 4. 

Finally, we present a dynamical phase diagram in Fig. 5, where we plot the driving frequency versus the ratio of the 
optical wavelength to the period of the substrate potential. The figure shows that, within a finite range of driving 
frequencies and wavelengths, the system as a whole will be running; otherwise, the particles will be locked at the 
bottom of the substrate potential, permanently. 

In conclusion, we would like to emphasize several points. First, a reasonable explanation for the photo-driven FK 
model seems to be that temporal and spatial modulations both result in the resonant absorption of optical energy. 
And, the dynamical competition between the inherent lengths of the moving object (the chain) and the supporting 
carrier (the isotropic surface) leads to the transformation of fed energy to directed motion. Second, the light-induced 
conformational changes of the photochromophore leads to the time-and-space dependence of the rest lengths of the 
bonds, which is the most important in this model. Third, according to our calculations, the steady state of the system 
is not affected by the initial conditions. The drift direction of the chain lies only on the sign in front of lu in the 
expression an,n+6 (t) = ao (1 -I- bcos{qXn,n+s T^t))- Namely, "— " and correspond to a motion in positive and 
negative a;-direction, respectively. Therefore, the motion of the system can be easily controlled. Fourth, the "springs" 
connecting the " particles" are light-sensitive molecules (photochromic molecules) and can expand or shrink with light 
reversibly JlSf ; adjusting the frequency or the wavelength of the lightwave may produce spectacular motions. We 
believe that the concept presented in this paper is simple and robust enough to be realized in actual experiments, 
and may demonstrate many interesting and typical nonlinear features of a system of interacting particles. However, 
we would like to further state that the connection between the model and photochromic chains stands at the level of 
speculation, being an interesting open question or idea, still to be proved. 
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Fig.l: Velocity of the centre of mass Vc of chain as a function of time t. For different parameters, the system turns 
quickly into a steady-going locked state ((a)-(c)) or a steady-going running state ((d),(e)) passing through a transitory 
relaxation. 

Fig.2: Average velocity of the centre of mass < Uc > of chain as a function of the amplitude b. 

Fig. 3: Average velocity of the centre of mass < Vc > of chain as a function of the driving frequency u. 

Fig. 4: Average velocity of the centre of mass < Wc > of chain as a function of A2/A1 that denotes the ratio of the 
optical wavelength to the period of the substrate potential. (Although the relation N/M = X2/X1 is not satisfied 
everywhere, the result is still credible with the large system size oi N = 1000.) 

Fig. 5: Dynamical phase diagram in the (w, A2/A1) plane for the photo-driven one-dimensional FK model. Here 
parameters are fixed at 6 = 0.3, rj = 0.1 and k ~ 1, respectively. And the particle-number N is chosen in 450 ~ 550 
to ensure that M waves (M is an integer, and M satisfy the relation: N/M = A2/A1) can be inserted into the 
commensurate background structure of A/'-particle system from first to last. 
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